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Abstract The Axelrod model is a spatial stochastic model for the dynamics of cultures which 
includes two important social factors: social influence, the tendency of individuals to become 
more similar when they interact, and homophily, the tendency of individuals to interact more 
frequently with individuals who are more similar. Each vertex of the interaction network is 
^— ■») ■ characterized by its culture, a vector of F cultural features that can each assumes q different 

' states. Pairs of neighbors interact at a rate proportional to the number of cultural features 

they have in common, which results in the interacting pair having one more cultural feature 
in common. In this article, we continue the analysis of the Axelrod model initiated by the 
first author by proving that the one-dimensional system fixates when F < cq where the slope 
satisfies the equation e"'^ = c. In addition, we show that the two-feature model with at least 
three states fixates. This last result is sharp since it is known from previous works that the 
one-dimensional two-feature two-state Axelrod model clusters. 

Oh 

, 1. Introduction 

^ . 

The Axelrod model is one of the most popular agent-based models of cultural dynamics. In addition 
to a spatial structure, which is modeled through a graph in which vertices represent individuals 
and edges potential dyadic interactions between two individuals, it includes two important social 
• factors: social influence and homophily. The former is the tendency of individuals to become more 

! similar when they interact, while the latter is the tendency of individuals to interact more frequently 

with individuals who are more similar. Note that the voter model [51 17] accounts for social influence 
since an interaction between two individuals results in a perfect agreement between them. The voter 
. model, however, excludes homophily. To also account for this factor, one needs to be able to define 

I a certain opinion or cultural distance between any two individuals through which the frequency of 

the interactions between the two individuals can be measured. In the model proposed by political 
scientist Robert Axelrod [T], each individual is characterized by her opinions on F different cultural 
^ I features, each of which assumes q possible states. Homophily is modeled by assuming that pairs 

of neighbors interact at a rate equal to the fraction of cultural features for which they agree, and 
social influence by assuming that, as a result of the interaction, one of the cultural features for 
which members of the interacting pair disagree (if any) is chosen uniformly at random, and the 
state of one of both individuals is set equal to the state of the other individual for this cultural 
feature. More formally, the Axelrod model on the one-dimensional lattice is the continuous-time 
Markov chain whose state space consists of all spatial configurations 
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that map the vertex set viewed as the set of all individuals into the set of cultures. To describe the 
dynamics of the Axelrod model, it is convenient to introduce 

F{x,y) ■■= p Yl = ?y(y,«)}, 

i=l 

where r]{x,i) refers to the ith coordinate of the vector r][x), which denotes the fraction of cultural 
features the two vertices x and y share. To describe the elementary transitions of the spatial 
configuration, we also introduce the operator (Jx,y,i defined on the set of configurations by 

|r/(y,i) if z = X and j = i 
[ax,y,iV)[z,j) '■= S , .u ■ for x,y G Z and z G |1,2, . . . 

I r][z,j) otherwise 

In other words, configuration dx^y^i i] is obtained from configuration ij by setting the ith feature of 
the individual at vertex x equal to the ith feature of the individual at vertex y and leaving the 
state of all the other features in the system unchanged. The dynamics of the Axelrod model is then 
described by the Markov generator L defined on the set of cylinder functions by 



2F 

\x-y\=l 4=1 



F{x,y) 
l-F{x,y) 



l{r]{x, i) 7^ r]{y, i)} [f{ax,y,i r?) - /(r?)]. 



The expression of the Markov generator indicates that the conditional rate at which the ith. feature 
of vertex x is set equal to the ith feature of vertex y given that these two vertices are nearest 
neighbors that disagree on their ith feature can be written as 
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F{x,y) 
1-F{x,y) 



F{l-F{x,y)) 2 



which, as required, is equal to the fraction of features both vertices have in common, which is the 
rate at which the vertices interact, times the reciprocal of the number of features for which both 
vertices disagree, which is the probability that any of these features is the one chosen for update, 
times the probability one half that vertex x rather than vertex y is chosen to be updated. Note 
that, when the number of features F = 1, the system is static, while when the number of states 
per feature q = 1 there is only one possible culture. Also, to avoid trivialities, we assume from now 
on that the two parameters of the system are strictly larger than one. 

The main question about the Axelrod model is whether the system fluctuates and evolves to a 
global consensus or gets trapped in a highly fragmented configuration. To define this dichotomy 
rigorously, we say that the system fluctuates whenever 

P {r]t{x, i) changes value at arbitrary large t) = 1 for all x G Z and i G {1, 2, . . . , F} (1) 

and fixates if there exists a configuration r]^ such that 

P {rit{x, i) = i]oo{x, i) eventually in t) = 1 for all x G Z and i G {1, 2, . . . , F}. (2) 

In other words, fixation means that the culture of each individual is only updated a finite number 
of times, so fiuctuation ([T|) and fixation ([2]) exclude each other. We define convergence to a global 
consensus mathematically as a clustering of the system, i.e.. 



lim P {r]t{x , i) = "qtiy , i)) = 1 for all x, ^/ G Z and i G {1, 2, . . . , F}. 



(3) 
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Note that whether the system fluctuates or fixates depends not only on the number of cultural 
features and the number of states per feature, but also on the initial distribution. Indeed, regardless 
of the parameters, the system starting from a configuration in which all the individuals agree for 
a given cultural feature while the states at the other cultural features are independent and occur 
with the same probability always fluctuates. On the other hand, regardless of the parameters, the 
system starting from a configuration in which all the even sites share the same culture and all the 
odd sites share another culture which is incompatible with the one at even sites always fixates. 
Also, we say that fluctuation/fixation occurs for a given pair of parameters if the one-dimensional 
system with these parameters fluctuates/fixates when starting from the distribution ttq in which 
the states of the cultural features within each vertex and among different vertices are independent 
and uniformly distributed. We also point out that neither fluctuation implies clustering nor fixation 
excludes clustering in general. Indeed, the voter model in dimensions larger than or equal to three 
for which coexistence occurs is an example of spin system that fluctuates but does not cluster while 
the biased voter model [21 [3] is an example of spin system that fixates and clusters. In spite of these 
counter-examples, we conjecture that fluctuation implies clustering and fixation excludes clustering 
for the one-dimensional Axelrod model starting from the distribution ttq. 

We now give a brief review of the previous results about the one-dimensional Axelrod model and 
state the new results proved in this article. Since two neighbors are more likely to interact as the 
number of cultural features increases and the number of states per feature decreases, one expects the 
phase transition between the fluctuation/clustering regime and the fixation/no clustering regime to 
be an increasing function in the F-q plane. The numerical simulations together with the mean-field 
approximation of [11] suggest that the system starting from ttq 

• exhibits consensus (clustering) when q < F and 

• gets trapped in a highly fragmented configuration (no clustering) when F < q. 

Looking now at analytical results, the first result in [^ states that the one-dimensional two-feature 
two-state Axelrod model clusters. The second result deals with the system on a large but finite 
interval, and indicates that, for a certain subset of the parameter region, the system gets trapped 
in a random configuration in which the expected number of cultural domains scales like the number 
of vertices. This strongly suggests fixation of the infinite system in this parameter region, which we 
prove in this paper. Shortly after, Lanchier and Schweinsberg [9] realized that the analysis of the 
Axelrod model can be greatly simplified using a coupling to translate problems about the model 
into problems about a certain system of random walks. To visualize this coupling, think of each 
spatial configuration as a q^-coloring of the set Z x {1, 2, . . . , F} and 

put a particle at {u,i) whenever rj{u — 1/2, i) ^ rj{u + 1/2, i) (4) 

for all u G Z + 1/2 and all cultural features i. We call u a blockade when it contains F particles, or 
equivalently when the two individuals on each side of u completely disagree. When the number of 
states per feature q = 2, Lanchier and Schweinsberg [9] proved that the construction (j4]) induces a 
system of annihilating symmetric random walks that has a certain site recurrence property, which is 
equivalent to fiuctuation of the Axelrod model, when starting from ttq. From this property, they also 
deduced extinction of the blockades and clustering, thus extending the first result of [8] to the model 
with two states per feature and any number of features. In contrast, the present paper deals with 
the fixation part of the conjecture and extends the second result of [^ by again using the random 
walk representation induced by (jl]). The first step is to prove that, for all values of the parameters. 
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the construction (j4|) induces a system of random walks in which cohisions result independently 
in either annihilation or coalescence with some specific probabilities. Coalescing events only occur 
when the number of states q > 2. This is then combined with large deviation estimates for the initial 
distribution of particles to obtain survival of the blockades when starting from ttq in the parameter 
region described in the second result of [8]. This not only implies fixation of the infinite system but 
also excludes clustering so the system gets trapped in a highly fragmented configuration. 



Theorem 1 - Assume that 

u;iq,F) := q(l-l] -F{l-- \ > 0. (5) 



Then, fixation ([2]) occurs and clustering ([3|) does not occur. 

Interestingly, though the second result in [8] relies on a coupling between the Axelrod model and 
a certain urn problem along with some combinatorial techniques that strongly differ from the 
techniques in our proof, both approaches lead to the same sufficient condition ([5]). The set of 
parameters described implicitly in condition ([5]) corresponds to the triangular set of crosses in the 
two diagrams of Figure [H which we obtained using a computer program. The picture suggests that 
this parameter region is (almost) equal to the set of parameters below a certain straight line going 
through the origin. To find the asymptotic slope, observe that ii F = cq then 



lim q ^uj(q,F) = lim (l ) — c(l ) 

r->oo g->oo ^ qj y qj 



e-" - c. 



q 

In other respects, if = c, then we have 



(cg-1) In(l-i) -ln(c) = (1 - cq) c 

OO °° C / 1 \ °° /I C 

^niq) ^n+liqj'^'^ n + 1 

n=l n=0 n=l ^ 
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from which we deduce that 



cq In - ^) > 111 (c) + In (^1 - and " ^) ^ ^ ~ g 

This proves that the condition in the theorem holds for F = cq and so all F < cq since oj is 
decreasing with respect to its second variable. In particular, fixation occurs whenever 

F < cq where c ~ 0.567 satisfies e~'^ = c. 

See Figure [T] for a picture of the straight line with equation F = cq. Finally, though a;(3, 2) = 
and therefore Theorem [1] does not imply fixation for the two-feature three-state Axelrod model, 
our approach can be improved to also obtain fixation in this case. 



Theorem 2 - The conclusion of Theorem^ holds whenever F = 2 and q = 3. 
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Figure 1. Phase diagram of the one-dimensional Axelrod model in the F-q plane. The diagram on the left-hand side 
is simply an enlargement of the diagram on the right-hand side that focuses on small parameters. The continuous 
straight line with equation _F = g is the transition curve conjectured in [11]. The set of crosses is the set of parameters 
for which the conjecture has been proved analytically: the vertical line of crosses on the left-hand side of the diagrams 
is the set of parameters for which fluctuation and clustering have been proved in [9] while the triangular set of crosses 
is the set of parameters such that uj{q,F) > for which fixation is proved in Theorem [1] The dashed line is the 
straight line with equation F = cq where the slope c is such that c = e~'^. 



Note that this fixation result is sharp since the first result in [8] gives fiuctuation and clustering of the 
two-feature two-state Axelrod model in one dimension. In particular, the two-feature model fixates 
if and only if the number of states per feature g > 2. To conclude, we note that, in contrast with the 
techniques introduced in [9] that heavily relies on the fact that the system starts from ttq, our proof 
of Theorem [T] easily extends to show that, starting from more general product measures, the one- 
dimensional system fixates under a certain assumption stronger than ([5]). However, the estimates of 
Lemmas [3] and [6l and consequently the condition for fixation, in this more general context become 
very messy while the proof does not bring any new interesting argument. Therefore, we focus for 
simplicity on the most natural initial distribution ttq. 

2. Coupling with annihilating-coalescing random walks 

As pointed out in [8], one key to understanding the Axelrod model is to keep track of the dis- 
agreements between neighbors rather than the actual set of opinions of each individual. When 
the number of states per feature q = 2, this results in a collection of non-independent systems of 
annihilating symmetric random walks. Lanchier and Schweinsberg [9] have recently studied these 
systems of random walks in detail and deduced from their analysis that the two-state Axelrod 
model clusters in one dimension. When the number of states per feature is larger than two, these 
systems are more complicated because each collision between two random walks can result in either 
both random walks annihilating or both random walks coalescing. In this section, we recall the 
connection between the Axelrod model and systems of symmetric random walks, and complete the 
construction given in [9] to also include the case g > 2 in which coalescing events take place. 
To begin with, we think of each edge of the graph as having F levels, and place a particle on 
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an edge at level i if and only if the two individuals that this edge connects disagree on their ith 
feature. More precisely, we define the process 

^t{u,i) := 1 {r]t{u- 1/2, i) j^r]t{u+ 1/2, i)} for all u G B := Z + 1/2 

and place a particle at site ti G B at level i whenever = 1. To describe this system, it is 

convenient to also introduce the process that keeps track of the number of particles per site, 

F 

Ct(n) := ^ (,t{u,i) for ah u G D, 

i=l 

and to call site u a j-site whenever it contains a total of j particles: Ct{u) = j. To understand 
the dynamics of these particles, the first key is to observe that, since each interaction between 
two individuals is equally likely to affect the culture of any of these two individuals, each particle 
moves one unit to the right or one unit to the left with equal probability one half. Because the 
rate at which two neighbors interact is proportional to the number of cultural features they have in 
common, a particle at {u, i) jumps at a rate that depends on the total number of particles located 
at site u, which induces systems of particles which are not independent. More precisely, since two 
adjacent vertices that disagree on exactly j of their features, and therefore are connected by an 
edge that contains a pile of j particles, interact at rate 1 — j/F, the fraction of features they share, 
conditional on the event that n is a j-site, each particle at site u jumps at rate 

r{j) := (l-^) J = ^--j forj/0, (6) 

which represents the rate at which both vertices interact times the probability that any of the j 
particles is the one selected to jump. Motivated by ([U]), the particles at site u are said to be active 
if the site has less than F particles, and frozen if the site has F particles, in which case we call 
u a blockade. To complete the construction of these systems of random walks, the last step is to 
understand the outcome of a collision between two particles. Assume that {u,i) and {u + are 
occupied at time t— and that the particle at (u, i) jumps one unit to the right at time t, an event 
that we call a collision and that we denote by 

(u, i) — > {u + at time t. 

This happens when the individual at x := n + 1/2 disagrees with her two nearest neighbors on her 
ith feature at time t— and imitates the ith. feature of her left neighbor at time t. This collision 
results in two possible outcomes. If the individuals at x and x + 1 agree on their ith feature just 
after the update, or equivalently the individuals at x — 1 and x + 1 agree on their ith feature just 
before the update, then {u + becomes empty so both particles annihilate, which we write 

{u, i) — ^ {u + 1, i) at time t. 

On the other hand, if the individuals at x and x + 1 still disagree on their ith feature after the 
update, then (u + l,i) is occupied at time t so both particles coalesce, which we write 

(u, i) — ^ (u + 1, i) at time t. 
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Figure 2. Coupling between the Axelrod model and annihilating-coalescing random walks. 



We refer to Figure [2] for an illustration of the coupling between the four- feature three-state Axelrod 
model and systems of annihilating-coalescing random walks. Each particle is represented by a cross 
and the three possible states by the colors black, grey and white. In our example, there are two 
jumps resulting in two collisions: an annihilating event then a coalescing event. We also refer the 
reader to Figure [3] for simulation pictures of the systems of random walks when F = 3. 

Lanchier and Schweinsberg [S] observed that, when q = 2, random walks can only annihilate, 
which was the key to proving clustering. This is due to the fact that, in a simplistic world where 
there are only two possible alternatives for each cultural feature, two individuals who disagree with 
a third one must agree. In our context, the individuals at x — 1 and x + 1 must agree just before 
the update when q = 2, which results in an annihilating event. In contrast, when the number of 
states per feature is larger, the three consecutive vertices may have three different views on their 
ith cultural feature, which results in a coalescing event. We point out that, since the system of 
random walks collects all the times at which pairs of neighbors interact, the knowledge of the 
initial configuration of the Axelrod model and the system of random walks up to time t allows 
to re-construct the Axelrod model up to time t regardless of the value of the parameters. There 
is however a crucial difference depending on the number of states. When q = 2, collisions always 
result in annihilating events so knowing the configuration of the Axelrod model is unimportant 
in determining the evolution of the random walks. In contrast, when q > 2, whether a collision 
results in a coalescing or an annihilating event depends on the configuration of the Axelrod model 
just before the time of the collision. The key to all our results is that, in spite of this dependency, 
collisions result independently in either an annihilating event or a coalescing event with some fixed 
probabilities. In particular, the outcome of a collision is independent of the past of the system of 
random walks though it is not independent of the past of the Axelrod model itself. 

To prove this result, we need to construct the one-dimensional process graphically from a perco- 
lation structure and then define active paths which basically keep track of the descendants of the 
initial opinions. First, we consider the following collections of independent Poisson processes and 
random variables: for each pair of vertex and feature (x, i) G Z x {1, 2, . . . , F}, 

• we let {Nx^i{t) : t > 0) be a rate one Poisson process, 

• we denote by Tx^i{n) its nth arrival time: Tx^i{n) := inf {t : Nx/i{t) = n}, 

• we let {Bx^i{n) : n > 1) be a collection of independent Bernoulli variables with 



P{BxM = +^) = P {Bx,i{n) = -1) = 1/2, 
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• and we let {Ux,i{n) : n > 1) be a collection of independent Uniform (0, 1). 

The Axelrod model is then constructed as follows. At time t = Tx^iin), we draw an arrow labeled i 
from vertex x to vertex y := x + Bj. i{n) to indicate that if 

< r{Ct^{u)) and CtAu)^0 where u = G B (7) 

then the individual at vertex y imitates the ith feature of the individual at vertex x. In particular, 
as indicated in ([6]), the rate at which the imitation occurs is equal to one half times the fraction 
of cultural features both vertices have in common divided by the number of features for which 
both vertices disagree, which indeed produces the local transition rates of the Axelrod model. The 
graphical representation defines a random graph structure, also called percolation structure, from 
which the process starting from any initial configuration can be constructed by induction based 
on an argument due to Harris [6]. Each arrow in this percolation structure is said to be active 
if condition ^ is satisfied. Note that whether an arrow is active or not depends on the initial 
configuration, and that the fact that an i-arrow from vertex x to vertex y at time t is active implies 
that the ith feature of y must be equal to the ith feature of x at time t. We say that there is an 
active i-path from (z, s) to {x, t) whenever there are sequences of times and vertices 

So = s < si < ■ ■ ■ < Sn+i = t and Xq = z, xi, . . . , Xn = X 

such that the following two conditions hold: 

1. For j = 1,2, ... ,n, there is an active i-arrow from Xj_i to Xj at time sj. 

2. For J = 0, 1, . . . , n, there is no active i-arrow that points at {xj} x (sj, Sj+i). 

We say that there is a generalized active path from (z, s) to (x, t) whenever 

3. For j = 1,2, . . . ,n, there is an active arrow from Xj^i to Xj at time Sj. 

Later, we will use the notations and to indicate the existence of an active i-path and a 
generalized active path, respectively. Conditions 1 and 2 above imply that 

for all (x, t) € Z X and all i, there is a unique 2 € Z such that {z, 0) ^ {x,t). 

Moreover, because of the definition of active arrows and simple induction, the ith cultural feature 
of vertex x at time t is equal to the initial value of the ith. cultural feature of z so we call vertex z 
the ancestor of vertex x at time t for the ith feature. In contrast, generalized active paths, which 
can be seen as concatenations of active i-paths for possibly different values of i, do not have such 
an interpretation, but the concept will be useful later to prove fixation. 

Lemma 3 - Conditional on the realization of the system of random walks until time t— and the 
event that [u, i) — > (n + 1, i) at time t, we have 

{u, i) — ^ {u + 1, i) at time t with probability {q — 1)^^ 

[u, i) — ^ (ti + 1, i) at time t with probability {q — 2) ■ [q — 1)^^ 
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(a) F = 3 and q = 2 (b) F = 3 and g = 3 (c) F = 3 and g = 4 

Figure 3. System of annihilating-coalescing random walks on the torus with 600 vertices. 



Proof. Let x := ii+ 1/2 € Z. Due to one-dimensional nearest neighbor interactions active i-paths 
cannot cross each other, from which we deduce that 

as{x — l,i) < as{x,i) < as{x + l,i) for ah s > (8) 

where as{ ■ ,i) denotes the ancestor at time s for the ith feature, i.e., 

{as{y, i),0) {y, s) for y E {x — l,x,x + 1} and ah s > 0. 

Moreover, conditional on the event of a cohision (n, i) — > {u + l,i) at time t, there is a particle at 
{u,i) and a particle at {u + l,i) at time t— therefore 

r]Q{at-{x±l,i)) = r]t-{x±l,i) / r]t-{x,i) = ??o(at-(x, i)) (9) 

From ([8|) and ([9]), we deduce that, conditional on {u,i) — > {u + at time t, 

as{x — l,i) < as{x,i) < as{x + l,i) for all s < t. 

In other respects, we have 

(ti, i) — ^ {u + at time t 

if and only if (n, i) — > [u + 1, i) at time t and i]t-{x — 1, i) = rjt~{x + 1, i) 

if and only if (n, i) — > {u + 1, i) at time t and r]o{at-{x — 1, i)) = rjQ{at^{x + 1, i)). 

In particular, the outcome - either an annihilating event or a coalescing event - of a collision at 
time t is independent of the realization of the system of random walks up to time t—. Moreover, 
since the initial states are independent and uniformly distributed, the conditional probability of an 
annihilating event is equal to the conditional probability 

P {X = Z \ X Z ^Y) (10) 

where X, Y, Z are independent uniform random variables over {1, 2, . . . , q}. By conditioning on the 
possible values of Y , we obtain that is equal to 

<? 1 
^ P{X = Z\Xj^jandZj^j)P{Y = j) = J2i{Q-'^)Q)-' = (^-l)"'- 
j=i j=i 
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Figure 4. Picture related to the proof of Lemma [l] Dashed lines represent active i-paths for some i whereas the 
continuous thick hue on the right-hand side is a generalized active path as defined in Section [5] 

Finally, since each collision results in either an annihilating event or a coalescing event, the condi- 
tional probability of a coalescing event directly follows. This completes the proof. □ 



3. Sufficient condition for fixation 

The main objective of this section is to extend a result of [1] to the Axelrod model, and obtain a 
sufficient condition for fixation which is based on certain properties of the active i-paths. 

Lemma 4 - For all (z, i) G Z x {1, 2, . . . , F}, let 

T{z,i) := inf {i : (z,0) (0,t)}. 
Then, the Axelrod model fixates whenever 



lim P {T{z, i) < oo for some z < —N and some i = 1,2, . . . , F) = 0. 



(11) 



Proof. Extending an idea of Bramson and Griffeath [3] and generalizing the technique in jlOj . we 
set Tj := for every cultural feature i and define recursively the sequence of stopping times 



inf {t > rij_i : rjt{0,i) / r?r,,,_i (0, «)} for j > 1. 



In other words, the stopping time Tij is the jth. time the individual at the origin changes the state 
of her ith. cultural feature. Also, for each cultural feature i, we define the random variables 



the ancestor of vertex at time r^j for the ith feature 



as well as the collection of events 

Bi := {rjj < oo for all j} and Gi^N 



{\aij \ < N for ah j}. 



See the left-hand side of Figure [H for a schematic illustration of the stopping times Tj j and the 
corresponding vertices ajj. The assumption (|lip together with reflection symmetry implies that, 
for each cultural feature i, the event Gi^N occurs almost surely for some A^. It follows that 

F F F F 

^(U^O ^ E ^(^') = E p{Br^([jGi^M)) = jZ p{[j{B,r^Gi,M)). 
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Since the event that the individual at the origin changes her culture infinitely often is also the event 
that at least one of the events occurs, in view of the previous inequality, in order to establish 
fixation, it suffices to prove that 

PiBiDGi^N) = for ah i G {1,2,... ,F} and ah iV > 1. (12) 

Our proof of (|12p relies on some symmetry properties of the Axelrod model that do not hold for 
the cyclic particle systems considered in [1]. First, we let 

It{x,i) := {z € Z : (x, i) is the ancestor of {z,i) at time t} 

be the set of descendants of {x,i) at time t, and denote by Mt{x,i) its cardinality. Since each 
interaction between two individuals is equally likely to affect the culture of each of these two 
individuals, the number of descendants of any given site is a martingale whose expected value is 
constantly equal to one. In particular, the martingale convergence theorem implies that 

lim Mt{x,i) = Moo{x,i) with probability 1 where £' |Moo(2;, i)| < oo 

therefore, for almost all realizations of the process, the number of descendants of (x, i) converges 
to a finite value. Since in addition the number of descendants is an integer-valued process, 

(j{x,i) := inf {t > : Mf (x, i) = Moo(x, i)} < oo with probability 1. 

Using that simultaneous updates occur with probability zero, we deduce that the set of descendants 
inherits the properties of its cardinality in the sense that, with probability one, 

lim It{x,i) = Ioo{x,i) and p{x,i) := ini {t > : It{x,i) = Ioo{x,i)} < oo (13) 

where, due to one-dimensional nearest neighbor interactions, Ioo{x,i) is a random interval which 
is almost surely finite. To conclude, we simply observe that, conditional on Gj.Ar, the last time the 
individual at the origin changes the state of her ith. cultural feature is at most equal to the largest 
of the stopping times p{x, i) for x € {—N, N) from which it follows that 

P {Bi n Gi^N) = P {p{x, i) = oo for some - N < x < N) = 

according to (jl3p . This proves (jl2p and therefore the lemma. □ 

4. Proof of Theorem [1] 

In view of Lemma HI in order to prove fixation, it suffices to show that the probability of the event 
in equation ([TT]) . that we denote by H^, tends to zero as — > oo. The first step is to extend 
the construction proposed by Bramson and Griffeath [4] to the Axelrod model, the main difficulty 
being that two active paths at different levels can cross each other. Let r be the first time an active 
i-path for some i = 1,2, ... ,F that originates from ( — oo, —N) hits the origin, and observe that 

T = inf {T{z, i) : z e (-oo, -N) and i = 1, 2, . . . , F} 

from which it follows that 

Hn := {r(2;, i) < oo for some (z, i) G (-00, -A) X {1, 2, . . . , F}} = {r < oo}. 
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Denote by z* < —N the initial position of this active path. Also, we set 

z- := mm{z eZ: {z,0) {0,t)} < z* < -N 

z+ := max{z € Z : (z,0) ~^ (Oi fo^^ some a < t} > 

and define / = We point out that z- < z* in general since vertex z* is defined from the set 

of active i-paths whereas vertex z^ is defined from generalized active paths that are concatenations 
of active i-paths with different values of i. See the right-hand side of Figure H] for an illustration 
where the two vertices are different. Now, note that each blockade which is initially in the interval 
/ must have been destroyed, i.e., turned into a set of F — 1 active particles through the annihilation 
of one of the particles that constitute the blockade, by time r. Moreover, active particles initially 
outside the interval / cannot jump inside the space-time region delimited by the two generalized 
active paths implicitly defined in (jl4|) . Indeed, assuming that such particles exist would contradict 
either the minimality of z_ or the maximality of z+. In particular, on the event if at, all the blockades 
initially in / must have been destroyed before time r by either active particles initially in / or active 
particles that result from these blockade destructions. To estimate the probability of this last event, 
we first give a weight of —1 to each particle initially active by setting 

(/)(u) := — Co(^) = —i whenever Co(^) = i F. 

To define 4>{u) when u is initially occupied by a blockade, we observe that by Lemma [3] the number 
of collisions required to break a blockade is geometric with mean q — 1. Moreover, each blockade 
destruction results in a total of F — 1 active particles. Therefore, we set 

(p{u) := ip{u) — [F — 1) whenever Co(^) = F 

where ip{u) are independent geometric random variables with mean q — 1. The fact that H]\{ occurs 
only if all the blockades initially in / are destroyed by active particles initially in / or active particles 
resulting from these blockade destructions, can then be written as 

if TV C < ^ 0(n) < I C < ^ (piu) < for some I < -N and some r > L (15) 
^ uei ^ ^ u=i ^ 



To understand the first inclusion, simply observe that the sum of the </>(tt) is equal to the number 
of collisions required to break all the blockades minus the total number of active particles initially 
in the interval I or created from the destruction of blockades initially in I. Since the number of 
collisions is bounded by the number of such active particles, all the blockades initially in I can only 
be destroyed if the number of such active particles exceeds the number of collisions required, which 
gives the first inclusion. The second inclusion simply follows from the fact that 

(-A^,0) C 2+) = I since z_ < -N and z+ > 0. 

The expression oiuj{q, F) can be understood heuristically as follows: since the (p{u) are independent, 
one expects that fixation occurs if E (j){u) > 0. But 

F-l 

EcPiu) = {Ei;{u)-{F-l))P{Co{u) = F) - J] iP(Co(u) = i) 

= {Ei;iu) + l)P{Co{u)=F) - iPiCo{u)=i) = qPiCo{u)=F)-ECo{u) 

i=0 
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which, since Coi^) = Binomial (F, 1 — l/q), is precisely equal to oj{q,F). To deduce rigorously 
fixation from the positiveness of the expected value, which is done in the next two lemmas, we now 
prove large deviation estimates for Hj\i. The first of these two lemmas will be used in the proof of 
the second one to show that the total number of collisions required to break all the blockades in a 
large interval does not deviate too much from its expected value. 

Lemma 5 - Let Xi, X2, ■ ■ ■ be an infinite sequence of independent geometric random variables 
witli the same parameter p. Tlien, for all e > 0, there exists 71 > such that 

P{Xi + X2 + --- + Xk <{l/p- e)K) < exp (-71K) for aU K sufficiently large. 

Proof. Let Z„ = Binomial (n,p) for all n > 1. Since, in a sequence of independent Bernoulli trials 
with success probability p, the event that the i^Tth success occurs at step n is included in the event 
that K successes occur in the first n steps, we have 

PiXi+X2 + ---+XK = n) < P{Zn = K). 

Letting M denote the integer part of {1/p — e)K, we deduce that 

M 

PiXi + X2 + --- + XK <{l/p-e)K) = P{Xi+X2 + --- + XK = n) 

M M "■^^ M 

< J2 P{Zn = K) < J2 P{Zn>K) < Y P{Zm>K) < MxP{Zm>K). 

n=K n=K n=K 

Since large deviation estimates for the binomial distribution imply that 

P{Zm>K) < P{ZM>{l-ep)~H'Ip) 

< exp(-72M) < exp(-72((l/p-e)is:- 1)) 

for a suitable constant 72 > 0, the result follows. □ 

Lemma 6 - Let In '■= {—N, 0) n ID and assume that uj{q, F) > 0. Then, 




for a suitable constant 73 > and all N sufRciently large. 
Proof. To begin with, we define 

Ni := card {u £ In : Co{u) = i} for z = 0, 1, . . . , F. 

Since the random variables Co(^)) "U G D, are independent, standard large deviation estimates for 
the binomial distribution imply that for all e > there exists 74 > such that 

P{N,^{{fii-e)N,{fi, + e)N)) < expi-j^N) for alH = 0, 1, . . . , F (16) 
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where := P (X = i) with X = Binomial (-F, 1 — l/q). The expression for fii fohows from the 
fact that initiaUy each level of each site is independently occupied with probability 1 — l/q, which 
implies that the Co('it) are independent binomial random variables. Let Q be the event that 

im - e)N <Ni< {fii + e)N for ah i = 0, 1, . . . , F. 

Then, there exists a constant C > such that, on the event il., 

^ F-l F-1 F 

^ iiVi < X] ^('"* + ^) - il^i-Ff^F + Ce = E Co{u) - F fip + Ce. 

i=0 1=0 i=0 

In particular, letting K be the integer part of (/x^? — e)A^, we have 

4>{u)<0\n\ < ^(V(^x)-(F-l)) < (i?Co(^^)-i^/UF + Ce)iv') 

< p(Y1 ^(^) ^ {ECo{u)-fiF + {C-F + l)e)N\ 



Now, since uj{q, F) > 0, there exists e > small such that 

ECo{u)-fiF + iC-F+l)e = {q-l)fiF + ECo{u)-qfiF + {C-F+l)e 

= {q-l)fiF-u}{q,F) + iC-F+l)e 
< {q-l-e){f,F-e) 

from which we deduce, also using ()17p and Lemma O that 

P( <P{u) < n\ < P( ^('") ^ - 1 - ) < exp (-71 if) 
for all K sufficiently large. Combining (|16|) and (|18p . we obtain 

P ( cPiu) < j < exp (-71 (/Xi. - e) TV) + (F + 1) exp {--f^N) 

for all N sufficiently large. □ 

Using the inclusion in p5]) and Lemma [H we deduce 



TV 

< 



(17) 



(18) 



lim P (Hm) < lim P ( > <i>(u) < for some / < —N and some r > | 

f~>OD ' ' ~ N->OD - / 

^ u=l ^ 

E E ^(E<^M<o) < Jin^ E exp(-73(r-/)) = 0. 

1<~N r>0 ^ «=/ ' 1<~N r>0 



This, together with Lemma HI implies fixation whenever uj{q,F) > 0. 
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5. Fixation when F = 2 and q = 3 

To begin with, note that, when F = 2 and q' = 3, we have E (f){u) = a;(3, 2) = for the comparison 
function (t){u) defined in the previous section. In particular, to find a good enough upper bound for 
the probabihty of in the case F = 2 and g = 3, one needs to define a new comparison function 
that also takes into account additional events that promote fixation, such as collisions between 
active particles and blockade formations. Recall that in the comparison function of Section HI each 
particle which is initially active is assigned a weight of —1, which corresponds to the worst case 
scenario in which the active particle hits a blockade. However, each active particle can also hit 
another active particle or form a new blockade with another active particle. More precisely, there 
are four possible outcomes for each active particle. 

1. If the active particle hits a blockade, it is assigned a weight of —1. 

2. If the active particle coalesces with another active particle, then at most one collision with 
a blockade can result from this pair of particles so the pair is assigned a total weight of — 1, 
i.e., each particle of the pair is individually assigned a weight of —1/2. 

3. If the active particle annihilates with another active particle, then no collision with a blockade 
can result from this pair so each active particle that annihilates with another active particle 
is assigned a weight of 0. 

4. If the active particle forms a blockade with another active particle, then following the same 
approach as in the previous section the pair is assigned a total weight equal to —1 plus a 
geometric random variable with mean (? — 1. 

In view of cases 2-4 above, the weight of an active particle that either collides with another active 
particle or forms a blockade with another active particle is at least —1/2, therefore we define a new 
comparison function, again denoted by (j), as follows: 

^(n)-l if Uu) = 2 

if CoM = o 

— 1/2 if Coi'^) = 1 ^iid the active particle initially at u either 

collides with another active particle or forms 
a blockade with another active particle 

— 1 if Co(^) = 1 find the active particle initially at u collides 
with a blockade 

where the random variables ^(tt) are again independent geometric random variables with the same 
expected value q — 1 = 2. The value of when Co('u) 7^ 1 is the same as in the previous section 
whereas we distinguish between active particles that satisfy case 1 or cases 2-4 above. The same 
reasoning and construction as in Section H] again imply that 

Hn C < < for some / < —N and some r > > (19) 

u=l J 

for this new comparison function. To prove that the probability of the event on the right-hand side 
converges to zero as iV ^> 00, we follow the same strategy as for Lemma [6] but also find a lower 
bound for the probability that a particle initially active either collides with another active particle 
or forms a blockade with another active particle, which is done in the next lemma. 



0(n) := < 
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Lemma 7 - Assume that F = 2 and q > 3. Then, there exists 75 > such that 




< exp (— 75A^) for all N sufRciently large 



where In '■= {—N, 0) n B as in Lemma El 

Proof. The first step is to find a lower bound for the initial number of active particles that will 
either collide or form a blockade with another active particle. To do so, we introduce the following 
definition: an active particle initially at site n € D is said to be a good particle if 

(0(14) = ((^(v) = 1 where {u,v} = {2n - 1/2, 2n + 1/2} for some n £ Z. (20) 

In other words, we partition the lattice D into countably many pairs of adjacent sites, and call 
an active particle at time a good particle if the other site of the pair is initially occupied by an 
active particle as well. An active particle which is not good is called a bad particle. Since initially 
each level of each site is independently occupied with probability 1 — 1/q, the variables Co(''^) are 
independent binomial random variables, so for u,v as in (I20p we have 

P {{u, v} is occupied by a pair of good particles at time 0) = = P {X = 1)^ 

where X = Binomial (2, 1 — l/q). Similarly, we have 

P (u is occupied by a bad particle at time 0) = i^i = P {X = 1) x P {X 7^ 1) 
P {u is occupied by two particles at time 0) = 1^2 = P (^ = 2). 

Since in addition the events that nonoverlapping pairs of adjacent sites are initially occupied by two 
good particles, or one bad particle, or one blockade, or one bad particle and one blockade, or two 
blockades are independent, standard large deviation estimates for the binomial distribution imply 
that there exists a positive constant 76 > such that 

PiN,^i{i^i-e)N,{i^, + e)N)) < exp{-jeN) fori = 0,1, 2, (21) 

where Nq,Ni and denote respectively the initial number of good particles, the initial number 
of bad particles, and the initial number of blockades in the interval In. To estimate the probability 
that a pair of good particles collide or form a blockade, we first observe that, when there are only 
two features, the graphical representation of the Axelrod model simplifies as follows: For each pair 
of neighbors {x,y) S Z^, draw an arrow x ^ y at the times of a Poisson process with intensity 
one fourth, which is equal to half of the rate at which neighbors who agree on one cultural feature 
interact. If the two neighbors agree on exactly one cultural feature at the time of the interaction 
then the culture of the individual at vertex y becomes the same as the culture of the individual at 
vertex x. In this graphical representation, there are exactly six possible arrows that may affect the 
system of random walks at the pair of sites {u,u + l}cB, namely 

u-1/2 u + 1/2 u + 3/2 u + 1/2 

n + 1/2 u-1/2 u + 1/2 u + 3/2 (22) 

u-3/2 u-1/2 u + 5/2 u + 3/2. 
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The event that one of the two arrows in the first hne of (I22p appears before any of the four other 
ones occurs with probabihty two (arrows) over six (arrows) = 1/3, and on the intersection of this 
event and the event that there is initially a pair of good particles at {u,u + 1}, the two particles 
either collide or form a blockade. Moreover, the event that one of the two arrows in the first line 
appears first only depends on the realization of the graphical representation in 



In particular, parts of the graphical representation associated with non-adjacent pairs do not inter- 
sect which, by independence of the Poisson processes, implies that the events that the two arrows 
in the first line of (j22p appears before any of the other ones are independent for non-adjacent pairs. 
It follows that the initial number J of good particles in J^v that either collide or form a blockade 
is stochastically larger than a binomial random variable with Nvq/2 trials and success probability 
one third. Large deviation estimates for the binomial distribution then imply that 



(n-3/2,ii-h5/2) X [0,oo). 



P (J < (1/6 - e)(z/o - e)N \ Nq > [uq - e)N) < exp ( 



77iV) 



(23) 



for a suitable constant 77 > 0. Now, let be the event that 



(i/j - e)N < Ni< {ui + e)N for i = 0, 1, 2 



and J>(l/6 



e){uo-e)N 



and observe that there exists a constant C > such that, on the event Q 



(l/2)J + (iVo + iVi- J) = iVo + iVi-(l/2)J 

< (z^o + + 2e)N - (l/2)(l/6 - e)(z.o - 



{111^0/12 + i^i + Ce)N. 



In particular, letting K be the integer part of (1^2 



e)A^, we have 




(24) 



In other respects, recalling the definition of Mj for i = 0, 1,2, we have 



{q -2)v2-v^- 11 1/0/12 

= ((? - 2) P (X = 2) - P (X = 1) P (X 7^ 1) - (11/12) P(X = 1)2 
= (g - 2) P (X = 2) - P (X = 1) (1/12) P{X = 1)2 



which, recalling the definition of X, is equal to 



(25) 
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for all g > 3. In particular, there exists e > small such that 

11 vqIYI + z^i + + (C - 1) e 

= (g - 1) 1^2 - ((g -2)v2-vx- 11 + (C - 1) e 

= (g - 1) 1^2 - (z^2 + g - 1) e < (g - 1 - e)(z^2 - e)- 

Since E'ip{u) = g — 1, the previous estimate, (|24p and Lemma [5] imply that 

p(Y1 '^(^) - ^ ) - ^( Yl < - 1 - e)^ ) < exp (-71 iT) 

for all sufficiently large. Combining H^U^ . ([25]) and we obtain 

^ ( "^(^^ - ^) ^ «^^P (^2 - e) A^) + 3 exp (-TeiV) + exp (-77iV) 

for all N sufficiently large, which completes the proof. □ 
As in the previous section, (|19|) and Lemma [7] imply that 

hm P{Hm) < lim Yl E exp (-75 (r-/)) = 

/<-7V r>0 

which, together with Lemma HI implies fixation when F = 2 and g = 3. 
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